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I.  Abstract 

For  high  Mach  number  flows,  M  >  4,  piston  theory  has  been  used  to  calculate  the  pressures  on  the 
surfaces  of  a  vehicle.  In  a  two-dimensional  inviscid  flow,  a  perpendicular  column  of  fluid  stays  intact  as  it 
passes  over  a  solid  surface.  Thus,  the  pressure  at  the  surface  can  be  calculated  assuming  the  surface  were  a 
piston  moving  into  a  column  of  fluid.  In  this  work,  first-order  piston  theory  is  used  to  calculate  the  forces, 
moments,  and  stability  derivatives  for  longitudinal  motion  of  a  hypersonic  vehicle.  Piston  theory  predicts  a 
relationship  between  the  local  pressure  on  a  surface  and  the  normal  component  of  fluid  velocity  produced 
by  the  surface’s  motion.  The  advantage  of  piston  theory  over  other  techniques,  such  as  Prandtl-Meyer  flow, 
oblique  shock,  or  Newtonian  impact  theory,  is  that  unsteady  aerodynamic  effects  can  be  included  in  the 
model.  The  unsteady  effects,  considered  in  this  work,  include  perturbations  in  the  linear  velocities  and 
angular  rates,  due  to  rigid  body  motion.  This  provides  a  more  accurate  model  that  agrees  more  closely  with 
models  derived  using  computational  fluid  dynamics  or  those  derived  by  solving  Euler  equations.  Additionally, 
piston  theory  yields  an  analytical  model  for  the  longitudinal  motion  of  the  vehicle,  thus  allowing  design  trade 
studies  to  be  performed  while  still  providing  insight  into  the  physics  of  the  problem. 

II.  Introduction 

In  the  1980’s,  the  National  Aerospace  Plane  (NASP)  program  commenced,  with  its  goal  being  a  feasibility 
study  for  a  single-stage  to  orbit  (SSTO)  vehicle,  which  was  reusable  and  could  take  off  and  land  horizontally. 
The  NASP  was  to  be  powered  by  a  supersonic  combustion  ramjet  (scramjet)  engine.  Although  this  program 
was  cancelled  in  the  1990’s,  a  great  deal  of  knowledge  was  gained  and  it  spawned  future  programs,  including 
the  hypersonic  systems  technology  program  (HySTP),  initiated  in  late  1994,  and  the  NASA  X-43A.  The 
HySTP’s  goal  was  to  transfer  the  accomplishments  of  the  NASP  program  to  a  technology  demonstration 
program.  This  program  was  cancelled  in  early  1995.  The  NASA  X-43A  set  new  world  speed  records  in  2004, 
reaching  Mach  6.8  and  Mach  9.6  on  two  separate  occasions  with  a  scramjet  engine.  These  flights  were  the 
culmination  of  NASA’s  Hyper-X  program,  with  the  objective  being  to  explore  alternatives  to  rocket  power 
for  space  access  vehicles. 

With  renewed  interest  in  space  operations  worldwide,  there  is  an  interest  in  hypersonic  aerodynamics 
research.  The  scramjet  engine  will  likely  play  a  major  role  in  future  hypersonic  vehicles.  Unlike  a  conventional 
turbojet  engine,  a  scramjet  engine  does  not  use  high  speed  turbomachinery  to  compress  the  air  before  it 
reaches  the  combustor.  Instead,  it  relies  upon  the  rise  in  pressure  across  oblique  shock  waves  located  in  front 
of  the  inlet.  Furthermore,  the  flow  through  the  entire  engine  is  supersonic  in  contrast  to  a  ramjet,  where  the 
flow  speeds  are  subsonic  through  the  combustor.  On  configurations  like  the  NASP  and  X-43A,  the  underside 
of  the  airframe  must  function  as  the  air  inlet  mechanism  and  the  exhaust  nozzle.  Therefore,  integration  of 
the  airframe  and  engine  are  critical  to  success  of  a  scramjet  powered  vehicle. 

Scramjets  could  be  used  as  part  of  a  multi-stage  launch  vehicle  that  would  include  multiple  propulsion 
systems  to  perform  a  mission.  The  factor  driving  research  towards  scramjets  and  away  from  rockets  is  cost; 
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scramjets  would  substantially  lower  costs  because  it  is  an  airbreathing  engine.  Airbreathing  engines  don’t 
require  oxidizer  to  be  carried  by  the  vehicle,  hence  increasing  the  payload  and  reducing  the  quantity  of  fuel 
carried. 

Unsteady  aerodynamics  is  another  key  technology  in  the  development  and  optimization  of  future  hyper¬ 
sonic  vehicles.  The  combined  effects  of  a  slender  flexible  vehicle  travelling  at  high  speeds  and  subjected  to 
large  forces  may  lead  to  significant  unsteady  aerodynamic  effects.  Hence,  understanding  the  concepts  and 
consequences  of  time-dependent  aerodynamic  flows  is  critical  to  the  success  of  this  type  of  vehicle. 

Piston  theory  is  a  technique  which  has  been  used  for  years  to  determine  the  pressure  distributions  on 
an  airfoil/ vehicle,  when  the  Mach  number  is  sufficiently  high.  Lighthill1  discussed  the  application  of  piston 
theory  on  oscillating  airfoils  some  50  years  ago.  Ashley  and  Zartarian2  discuss  piston  theory  while  providing 
a  number  of  examples  of  the  application  of  piston  theory  to  specific  problems.  More  recently,  Tarpley3 
discussed  the  computation  of  stability  derivatives  for  a  caret-wing  waverider  using  piston  theory,  which 
requires  the  analysis  of  unsteady  flow  over  the  vehicle.3  Piston  theory  allows  the  inclusion  of  unsteady 
aerodynamic  effects  in  the  model  and  a  closed  form  solution  can  be  found  for  these  unsteady  effects. 

In  this  work,  piston  theory  is  applied  to  a  hypothetical  2-dimensional  hypersonic  vehicle  powered  with 
a  scramjet.  Section  III  describes  the  vehicle  analyzed  in  this  work,  while  the  steady  forces  on  the  vehicle’s 
surfaces  are  calculated  in  Section  IV.  The  afterbody  effects  are  included  in  Section  V,  where  the  pressure 
distribution  and  force  due  to  exhaust  plume  are  evaluated.  The  control  surface,  which  is  a  single  elevator, 
is  examined  in  Section  VI,  analysis  of  the  flow  regions  is  performed  in  Section  VII,  total  body  forces  and 
moments  are  derived  in  Section  VIII,  while  stability  derivatives  are  calculated  in  Section  IX.  The  engine 
model  is  developed  in  Section  X  and  some  results  from  an  open-loop  simulation  are  provided  in  Section  XIII. 

III.  HSV  Model 

Figure  1  shows  the  hypersonic  vehicle  considered  in  this  work.4  The  vehicle  consists  of  4  surfaces: 
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Figure  1.  Hypersonic  Vehicle. 


an  upper  surface  (defined  by  points  cf)  and  three  lower  surfaces  (defined  by  points  cd,  gh,  and  ef).  All 
pertinent  lengths  and  dimensions  are  in  units  of  feet  and  degrees,  respectively.  The  total  length  of  the 
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vehicle  is  L  =  100 ft  and  the  notation  for  lengths  is  Lf  —  length  of  the  forebody,  Ln  =  length  of  the  engine 
nacelle,  La  =  length  of  the  aftbody,  Le  is  the  length  of  the  elevator,  Xf  is  the  distance  from  the  C.G.  to  the 
front  of  the  vehicle,  xa  is  the  distance  from  the  C.G.  to  the  rear  of  the  vehicle,  xcs  and  zcs  are  the  distances 
from  the  C.G.  to  the  midpoint  of  the  elevator  in  the  x  and  z  directions,  respectively,  and  hi  is  the  engine 
height.  The  vehicle  lengths  are 

L  =  100  ft 
Lf  =  47  ft 
La  =  33  ft 
Ln  =  20  ft 
Le  =  17  ft 
Xf  =  55  ft 
Xa  =  45  ft 
Xcs  =  30  ft 
zcs  =  3.5  ft 
hi  =  3.5  ft 

The  vehicle  angles  are 

ri,u  =  3° 

7-1, L  =  6°  (2) 

r2  =  14.41° 

Additionally,  the  vehicle  mass  and  moment  of  inertia  are 

Mass  =  300  slug  ,  * 

Jyy  =  500, 000  slug  —  ft 2 

and  the  mean  aerodynamic  chord  (c)  and  planform  area  (. S )  are  defined  as 


c  =  L 
S  =  L 2 


(4) 


The  goal  is  to  apply  piston  theory  to  this  vehicle  to  determine  the  pressure  distribution  on  the  surfaces 
of  the  vehicle,  which,  in  turn,  can  be  used  to  evaluate  the  forces  and  moments.  The  pressure  on  the  face  of 
a  piston  moving  into  a  column  of  perfect  gas  is2 


7-1  Vn\ 

2  a0 o  J 


27 

7-1 


(5) 


where  the  subscript  ”oo”  refers  to  the  steady  flow  conditions  past  the  surface,  Vn  is  the  velocity  of  the  surface 
normal  to  the  steady  flow,  is  the  freestream  speed  of  sound,  and  P  is  the  surface  pressure.  Taking  the 
binomial  expansion  of  Eq.  5  produces 


P_  =  x  +  2 7  7  -  1  Vn  =  1  +  7  Vn 

Poo  7  1  2  Q/0C  ^oo 


(6) 


Multiplying  through  by  P^  and  using  the  perfect  gas  law  (P  =  pRT )  and  the  definition  of  the  speed  of 
sound  (a2  =  7 RT)  yields  the  basic  result  from  first-order  linear  piston  theory 


P  —  Poo  T  Poo^oo  bn 


(7) 


where  7  is  the  ratio  of  specific  heats  and  R  is  the  gas  constant.  The  infinitesimal  force  due  to  the  pressure  is 

dF  =  -PdAn  (8) 

where  dA  is  a  surface  element  and  n  is  the  outward  pointing  normal.  Substituting  Eq.  7  into  Eq.  8  yields 

dF  =  (-Poo  -  PoodooVn)  dAn  (9) 
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The  normal  velocity  can  be  computed  by  taking  the  dot  product  of  the  flow  velocity  over  a  surface  and  the 
outward  pointing  normal  for  that  surface.  Hence,  Eq.  9  becomes 

dF  =  (-Poo  -  Poodoc  [v  ■  n])  dAn  (10) 

Equation  10  is  the  basic  result  upon  which  this  work  is  based.  From  this  equation,  it  is  seen  that  in  order 
to  compute  the  forces  acting  on  a  surface,  one  must  determine  the  properties  of  the  flow  past  the  surface 
(properties  behind  a  shock,  expansion  fan,  or  freestream),  the  velocity  of  the  surface  relative  to  the  airstream, 
V,  the  outward  pointing  surface  normal,  n,  and  the  surface  element,  dA.  The  work  that  follows  will  develop 
these  quantities  for  the  upper  surface  and  the  lower  surfaces. 

IV.  Vehicle  Surface  Pressure  Distributions  and  Forces 

To  compute  the  forces,  moments,  and  stability  derivatives,  consider  small  perturbations,  from  a  steady 
flight  condition  at  M^,  in  the  velocities  u  and  w  and  the  rate  q.  On  the  upper  surface,  the  surface  is  modelled 
as  a  piston  moving  into  a  column  of  fluid  that  has  the  properties  of  the  flow  behind  an  oblique  shock  wave,  an 
expansion  fan,  or  freestream  (flow  properties  are  determined  by  the  angle  of  attack).  Likewise,  on  the  lower 
surface,  the  surface  is  modelled  as  a  piston  moving  into  a  column  of  fluid  that  has  the  properties  of  the  fluid 
behind  the  oblique  shock.  An  oblique  shock  is  required  on  the  lower  surface  for  proper  engine  operation. 
Figure  2  shows  the  regions  of  interest. 


Figure  2.  Hypersonic  Vehicle,  Oblique  Shocks,  and  Pressure  Regions. 


To  begin  the  development,  first  consider  the  upper  surface.  The  velocity  of  a  point  on  the  upper  surface 
due  to  the  velocity  and  rate  perturbations  is 

Vc/  =  ( Vi  cos  T\  u  +  u)  i  +  ( Vi  sinrpc/  +  w)  k  +  cj  x  rc/  (11) 

where  z,  k  are  unit  vectors  in  the  x  and  z  body  axes,  respectively,  uj  is  the  angular  rate  vector,  a  is  the  angle 
of  attack,  and  V\  is  the  velocity  of  the  flow  on  the  upper  surface  (see  region  1  in  Fig.  2).  For  longitudinal 
motion  only,  gj  =  qj  where  j  is  a  unit  vector  in  the  y  body  axis  direction.  In  Eq.  11,  r cf  is  the  position 
vector  of  a  point  on  the  upper  surface  given  by 

r cf  =  rcfj  +  +rcfzk  xi  +  tanri?c/  (x-Xf)k 

—  Xa  <  X  <  Xf 
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According  to  Eq.  8,  a  normal  vector  to  the  upper  surface  is  also  needed.  The  upper  surface  outward  pointing 
normal  vector  is 

ncf  =  sinTi?tA  —  cos  T\jjk  (13) 

For  the  lower  surface  defined  by  the  points  c  and  d  in  Figure  1,  we  use  the  velocity  of  the  flow  after  the 
oblique  shock  to  obtain 

Vcd  =S  (V2  cos  ti5l  +  u)  i  +  (- V2  sin  r1}L  +  w)  k  +  uj  x  rcd  (14) 

while  for  the  surface  defined  by  points  g  and  h 

V gh  =  (V3  +  u)  i  +  wk  +  cv  x  rgh  (15) 

where  r cd  and  rgh  are  position  vectors  of  a  point  on  the  lower  surface  given  by 

rcd  =  rcdJ  +  rcdk  =  xi-  tanri?L  (x-xf)k  ^ 

Xf  —  Lf<x<Xf 

r gh  =  rghJ  +  rghzk  =  xi  +  (Lf  tan  74, L  +  hi)  k  n  ^ 

(xf  —  Lf)  —  Ln  <  x  <  Xf  —  Lf 

The  normal  vectors  for  the  lower  surfaces  are 

ncd  =  sinri^i  +  cosri  ,Lk 
ngh  =  lk 

Performing  the  cross  products  required  by  Eqs.  11,  14,  and  15  gives 

c o  x  r cf  =  qtsmri^u  (x  —  Xf)i  —  qxk 

uj  xvcd  —  —  gtanri?L  (x  —  Xf)i  —  qxk 
oj  x  r  gh  —  q[Lf  tanri^  +  hi)  i  —  qxk 
According  to  Eq.  10,  the  pressures  on  the  surfaces  of  interest  are 

Pcf  —  Pi  T  Pi^i  (' Vc/  •  n cf) 

Pcd  =  P2  +  P2CL2  (Vcd  •  ncd) 

Pgh  =  P3  +  P3CL3  (Vgh  ■  n^) 

Substituting  the  results  of  Eq.  22  into  Eq.  8  gives 

dF cf  —  {  Pi  —  Picti  (Vc/  •  n c/)}  dAcfUcf 
dF cd  —  {  P2  P2&2  (V cd  '  ^cd)}  dAcdYlcd 
dFgh  =  {-P3  -  p3a3  (' Vgh  ■  n^)}  dAgh ngh 

Using  Eqs.  11,  14,  and  15  and  the  appropriate  normal  vectors  (Eqs.  13  and  18),  the  dot  products  in  Eq.  23 
become 

Vc/  •  n cf  =  [u  +  gtanri?[/  (x  —  Xf)]  sinri^u  —  [w  —  qx]  costi?[/ 

Vcd  ■  ncd  =  [u-  gtanri?L  (x  -  xf)\  sin  74, L  +  [w  -  qx]  cos  rliL  (24) 

Vgh  ■  ngh  =  w  -  qx 

Note  that  the  steady  terms  cancel  as  a  result  of  taking  the  dot  product.  Using  Eq.  24  in  Eq.  23  yields 

dFcf  =  {—Pi  —  Picti  {[u  +  qtSiUTi^u  (x  —  Xf)]  sinri?[/  —  [w  —  qx]  cos  r^u})  dAcfncf 

dFcd  =  {-P2  -  p2a2  {[u  -  qtanrliL  (x  -  xf)]  sinri,L  +  [w  -  qx]  cos rliL})  dAcdncd  (25) 

dFgh  =  {-P3  -  p3a3  {w  -  qx})  dAghngh 


(18) 

(19) 

(20) 
(21) 

(22) 


(23) 
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The  next  step  is  to  determine  the  upper  and  lower  surface  elements.  Note  that  the  vehicle  model  is  2- 
dimensional  with  unit  depth  into  the  page.  Hence,  the  upper  surface  element,  dAcf  can  be  written  as 

dAcf  =  dLcf{  1)  (26) 

where  dLcf  defines  a  length  of  interest  on  the  upper  surface  and  the  multiplying  factor  of  1  is  due  to  the 
vehicle’s  unit  depth.  The  surface  element  can  be  written  as 

dAcf  =  \J  dx 2  +  dz2(  1)  (27) 


From  Eq.  12, 


Using  Eq.  28  in  Eq.  27  yields 


z  =  tanri?[/  (x  —  Xf)  =>  dz  =  tsaiTi^dx 


dAcf  =  y/ dx 2  +  tan2  ri^jjdx2  (1)  =  dxyjl  +  tan2  ~ti~u{1)  =  dx  secri?t/(l)  =  secTi^udx 
Similarly,  the  surface  elements  for  the  lower  surfaces  become 


(28) 

(29) 


dAcd  =  secTi^dx  (30) 

dAgh  —  dx  (31) 

Now,  the  incremental  force  on  the  upper  surface  (see  first  entry  in  Eq.  25)  becomes 

dFcf  =  (—Pi  —  picti  {[u  +  qt&nTi,u  (%  ~  %f)]  sinri :u  —  [w  —  qx]  cosri^})  ncf  secri^dx  (32) 

Using  similar  analysis,  the  incremental  forces  on  the  lower  surfaces  become 


dFcd  =  (-P2  -  P2d2  {[ u  -  gtanrpL  (x  -  Xf)\  sinrpz,  +  [w  -  qx]  cos r1?L})  ncdsec r\^dx  (33) 

dFgh  =  (-P3  -  p3a3  {w  -  qx})  nghdx  (34) 

It  will  be  useful  to  separate  the  steady  forces  from  the  unsteady  forces  at  this  point.  In  Section  IX,  we 
will  model  the  unsteady  effects  using  a  stability  derivative  approach.  To  compute  the  steady  forces,  the 
components  of  the  incremental  forces  related  to  the  steady  flow  are  integrated  over  the  surface  of  the  vehicle. 
For  the  upper  surface,  this  becomes 


■  cf 


' steady 

while,  for  the  lower  surfaces,  the  forces  are 

rxf 

=  (' d^cd)  steady 

Jxf-Lf 


fxf  fxf  r  „  „-i 

/  (dFcf) steady  =  /  -Pi  sin rljC7«  -  cos r1>Crfe  seen tUdx 

J—  Xn  J—  Xn  L  J 


rxf 

J  Xf~. 


xf~Lf 


sinri^i  +  cos  T\^k 


sec  r\  Ldx 


■  gh 


I, 


xf~Lf 


(dFgh) 


rxf~Lf 
J(xf-Lf 


steady 


—P3kdx 


(35) 


(36) 


(37) 


(xf  Lf)  Ln  Pxf  Lf)  Ln 

For  the  consideration  of  aeroelastic  effects  in  future  work,  the  upper  surface  force  is  separated  into  forebody 
and  afterbody  components.  Thus,  Eq.  35  becomes 


F  c/  =  /  -A 

~  Pcfa  +  Fcff 

Performing  the  integrations  yields 


sin  r\jji  —  cos  r\^jk  sec  T\fjdx  + 


rxf 

Jo 


-Pi 


sinri  xji  —  cosri  uk 


sec  n  udx 


(38) 


F cfa  =*  -PiXa  secret/ 
Fcff  =  —P\Xf  secri?i/ 


sinri5[/i  —  cosTi^uk 
sin  r\  jji  —  cos  r\\jk 


~  XCfa  i  +  Zcfa  k 
—  X-cffi  T  Zcff  k 


(39) 
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where  Xcfa ,  Zcfa  are  the  components  of  the  aftbody  upper  surface  force  in  the  x  and  z  directions,  respectively, 
and  Xcf  ,  Zcff  are  the  components  of  the  forebody  upper  surface  force  in  the  x  and  z  directions.  For  the 
lower  surfaces,  integration  of  Eqs.  36  and  37  produce 


F cd  =  -PzLf  sec  ti:L 


sinri  +  costi  lP 


—  Xcd ^  T  Zcdk 


F  gh  PsLnk  Xghi  +  Zghk 


(40) 


where  X *  and  Z *  are  the  components  of  the  lower  surface  forces.  Equations  39  and  40  give  the  rigid  body 
steady  forces  on  this  vehicle  due  to  pressure  distributions  on  the  upper  and  lower  surfaces. 


V.  Afterbody 


The  flow  on  the  afterbody  of  the  vehicle  is  bounded  by  the  vehicle  surface  (surface  ef  in  Fig.  1)  and  a 
shear  layer  between  the  freestream  atmosphere  and  the  exhaust  gas  of  the  engine.  Hence,  external  nozzle 
analysis  must  be  performed  to  determine  the  pressures  on  surface  ef.  According  to  Chavez,5  the  pressure 
distribution  along  the  nozzle  surface  can  be  approximated  by 

Pef(s)  » - - 

1  +  — £— 

cos{Tl,U+T2 

where  Pef  is  the  pressure  on  the  afterbody,  Pe  is  the  pressure  at  the  engine  exit,  P ^  is  the  freestream 
pressure,  and  s  is  the  distance  from  the  lower  apex  (point  e)  to  the  point  of  interest  along  the  vehicle’s 
afterbody  surface  (see  Fig.  1).  The  force  produced  by  the  external  nozzle  can  be  calculated  by  integrating 
Eq.  41  over  the  rear  ramp  of  the  vehicle: 


Fef  = 


I; 


s(Tl,t/+T2) 


1  + 


cos(Tl,C/  +  r2) 


(ft-1) 


ds 


Performing  the  integration  and  simplifying  yields 


Fef  = 


l-n  ft  l’\ 


COS  (nif/  +  72 )  ( Pe  Poo) 


In- 


(42) 


(43) 


Equation  43  provides  the  magnitude  of  the  force  due  to  the  external  nozzle.  Its  direction  is  perpendicular 
to  the  rear  ramp.  Hence,  the  vector  force  due  to  the  external  nozzle  is 


■ef 


LaPePoi 


COS  (71^+72)  (. Pe  Poo  ) 


ln- 


sin  (ntu  +  7-2)  i  —  cos  (r1:U  +  r2)k  =  Xefi  +  Zefk  (44) 


where  Xef  and  Zef  are  the  axial  and  normal  force  components  of  the  external  nozzle  force. 

For  use  in  calculating  stability  derivatives,  it  is  necessary  to  determine  the  force  on  the  rear  ramp  due  to 
perturbations  in  the  velocities  u  and  w  and  the  rate  q.  The  differential  force  on  the  rear  ramp  is 


dP ef  —  (  Pef  Pef  f  Wef  '  **e/])  dA.efXlef 


(45) 


where  Pef  is  the  afterbody  pressure  distribution  given  by  Eq.  41,  pef,aef  are  the  density  and  speed  of  sound 
on  the  rear  ramp,  Vej  is  the  flow  velocity  on  lower  surface  ef,  n ef  is  the  normal  vector  to  lower  surface  ef, 
and  dAef  is  the  surface  element  for  surface  ef.  The  position  vector,  normal  vector,  and  surface  element  for 
this  surface  are  given  by: 


r ef  —  xi  +  [tan  (t^u  +  t-2)  (x  +  xa)  -  Ltanr^c/]  k  =  xi  +  refk 
Xa  T  x  ^  La  xa 


(46) 


n ef  =  -  sin  (ti?C7  +  r2)  i  +  cos  (r1:U  +  r2)  k  (47) 

dAef  =  sec  (ri}c/  +  r2)  dx  (48) 

Using  Eqs.  46-48  in  Eq.  45  yields 

dFef  =  {-Pef  -  PefCLef  {-  {u  +  qrefz)  sin  +  t2)  +  (w  -  qx)  cos  (t^u  +  7-2)})  nefdAef  (49) 

Integration  of  the  steady  component  of  Eq.  49  yields  the  same  result  as  in  Eq.  44.  The  unsteady  components 
will  be  utilized  later  in  this  work  for  the  computation  of  stability  derivatives. 
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VI.  Control  Surfaces 


The  control  surface  is  an  elevator  located  near  the  tail  of  the  vehicle  as  shown  in  Fig.  1.  The  elevator  is 
modelled  as  a  flat  plate  hinged  at  its  midpoint  so  the  entire  surface  deflects.  The  length  of  the  elevator  is 
Le  =  17 ft.  Positive  8e  is  defined  as  trailing  edge  down.  Once  again,  the  velocity  of  the  flow  on  both  sides  of 
the  elevator  must  be  determined.  Proceeding  in  a  manner  similar  to  that  which  has  already  been  done,  it  is 
found  that 


and 


V eL  =  (V4  cos  5e  +  u)  i  +  (— V4  sin  5e  +  w)  k  +  gj  x  re 


Veu  =  (V5  cos  5e  +  u)  i  +  (— V5  sin  5e  +  w)  k  +  gj  x  re 


(50) 


(51) 


where  VeL  is  the  flow  velocity  on  the  underside  of  the  elevator,  Veu  is  the  flow  velocity  on  the  upper  surface 
of  the  elevator,  c 0  =  qj,  re  is  a  position  vector  from  the  vehicle  c.g.  to  an  arbitrary  point  on  the  elevator,  and 
V4,  V5  are  fluid  velocities  (freestream,  behind  oblique  shock,  or  behind  expansion  fan).  The  position  vector 
is  found  to  be 

ve  —  xi—  [ zcs  -f  tan  Se  (x  +  xcs)\  k 
~XCS  —  ^2  COS  <  x  <  —Xcs  +  cos  5e 

where  xcs  and  2;cs  are  the  x  and  z  positions  of  the  midpoint  of  the  elevator  referenced  to  the  c.g.  As  shown 
in  Fig.  1,  xcs  —  —30 ft  and  zcs  =  —3.5 ft.  For  this  control  surface,  outward  pointing  normal  vectors  for  both 
the  lower  and  upper  surfaces  are  needed.  These  normals  are  computed  as 


(52) 


neu  =  —  sin  5ei  —  cos  Sek 
neL  =  sin  Sei  +  cos  Sek 

To  evaluate  the  cross-product  in  Eq.  51,  we  use  Eq.  52  to  obtain 

wxr  e  —  —q  [ zcs  +  tan  Se  (x  —  xcsj]  i  —  qxk 

Then,  the  differential  forces  on  the  upper  and  lower  surfaces  of  the  elevator  become 

dFeu  —  [  7^5  P5^5  {Vef/  ■  neudAe 
dFeL  —  [  P4  P4.CL4.  {VeL  •  nejL|]  n eLdAe 


where 


Vet7  ■  neu  =  -  (u  -  q  {zcs  +  tan  5e  (x  -  xcs)})  sin  5e  -  (w  -  qx)  cos  Se 


Vei  ■  neL 


=  (u  —  q  {zcs  +  tan  Se  {x  —  xcs)})  sin  Se  +  (w  —  qx)  cos  Se 


and 


dAe  —  sec5edx(l)  —  sec  5edx 

Hence,  the  upper  and  lower  forces  on  the  elevator  can  be  computed  as 


-l 


-xcs-\-^f-  cos Se 


X cs  o  COS  Se 


-P5  -  p5a5  {Veu  .  neu  }] 


[~Pa  ~  PAaA  {V eL  .  nei}] 


•  sin  8  ei  —  cos  5Pk 


sin  8ei  +  cos  8ek 


sec  5Pdx 


:  8Pdx 


(53) 

(54) 

(55) 

(56) 

(57) 

(58) 

(59) 

(60) 


Using  the  steady  component  of  Eq.  56  in  Eq.  59  and  the  steady  component  of  Eq.  57  in  Eq.  60,  the  steady 
forces  on  the  elevator  become 


L 

-L 


cos  Se 

—xcs-\-^g  cos 8e 
—xcs  —  hf-  cos  8e 


-  sin  8 ei  —  cos  8ek  sec  8Pdx 


-Pa 


sin  8ei  +  cos  5ek 


:  8Pdx 


(61) 

(62) 
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Evaluating  these  integrals  yields 


Feu  =  —  P$Le  —  sin  5ei  —  cos  Sek  —Xeui-\-Zeuk 


Fer  =  —P^Le  |sin(5ei  +  cos5e^|  =  XeLi  +  ZeLk 


VII.  Flow  Analysis 


(63) 

(64) 


In  the  preceding  analysis,  the  properties  of  the  flow  on  the  upper,  lower,  and  control  surfaces  have  not 
been  defined.  In  this  section,  the  properties  of  the  flow  will  be  determined.  Specifically,  the  angles  of 
at  tack/ control  deflections  at  which  shock  waves  or  expansion  fans  are  created  will  be  delineated. 

By  examination  of  Fig.  1,  the  following  relationships  can  be  determined: 


if  a  =  rip  ->  freestream  ’•V\  =  V00,pi=p00,a1=aQ 
if  ol  >  — y  expansion  fan 

if  a  <  t\jj  shock(compression) 


(65) 


The  wedge  angles  associated  with  the  upper  surface,  for  calculation  of  flow  properties  behind  the  shock  or 
expansion  fan,  are  as  follows: 

if  a  >  Titu  — >•  0Uexpansion  =a-rhU 
if  a  <  T\tu  — >  0Ushock  —  -a  +  rltU 


(66) 


The  above  information  is  used  to  determine  the  flow  properties  for  the  upper  surface,  namely,  Vi,pi,  and 
a\.  For  lower  surface  cd,  the  relationships  become 


if  ol  —  ~rlyL  freestream  :  V2  =  Foo,  P2  =  Poo,  a2 
if  ol  >  —ti:l  -a  shock(compression) 
if  a  <  —ti:l  expansion f an 


(67) 


Lower  surface  cd  wedge  angles  are  as  follows: 


if  a  >  -TitL  0L 

shock  a  +  T\tL 


if  a  <  -Ti'L  >  9i 


'  expansi 


■  -OL  ~  TiiL 


(68) 


Physically,  for  the  scramjet  engine  to  work  properly,  an  oblique  shock  must  form  on  the  underside  of  the 
vehicle  to  increase  the  pressure  at  the  inlet.  This  effectively  places  a  lower  limit  on  the  angle  of  attack  and 
requires  that  the  angle  of  attack  be  such  that  an  oblique  shock  forms  on  the  underside  of  the  vehicle.  Hence, 
from  Eq.  67,  the  angle  of  attack  must  satisfy 


ol  >  Tl,l 


(69) 


In  other  words,  —  r\^  is  an  absolute  lower  limit  for  angle  of  attack,  with  the  engine  on;  however,  the  engine 
will  cease  to  function  before  this  lower  limit  is  reached.  With  this  limit  in  place,  a  bow  shock  will  form  on 
the  underside  of  the  vehicle.  For  lower  surface  gh,  it  is  necessary  to  calculate  the  angle  at  which  the  shock 
exactly  impinges  on  the  point  g  of  the  engine  nacelle.  This  angle,  denoted  by  Tbowshock ,  is 


Tbowshock  —  ol  T  tan 


-1 


Lf  tan  TpL  +  hi 


(70) 


Let  the  bow  shock  angle  be  denoted  by  (3.  If  (3  >  r^ow shock ,  then  the  shock  misses  the  point  g  and  the  flow 
properties  on  lower  surface  gh  are  computed  using  the  flow  properties  behind  the  oblique  shock  (F2,p25^2) 
as  the  initial  conditions.  Then,  an  expansion  fan  forms  at  point  g.  If  (3  <  Tbowshocki  then  either  the  shock  is 
on  the  lip  (point  g)  or  the  shock  is  inside  the  engine  inlet.  In  either  case,  freestream  properties  are  used  to 
compute  Vs ,  ps ,  a% .  The  following  steps  are  used  to  determine  the  flow  conditions  over  surface  gh. 


1.  Calculate  shock  angle,  /?,  from  surface  cd. 
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2.  If  (3  >  Ti)OWShocki  an  expansion  fan  forms  at  point  g  and  the  flow  properties  behind  the  oblique  shock 
wave  are  used  as  input  to  the  expansion  fan  flow  equations  with  a  wedge  angle  of  t\,l. 

3.  If  (3  <  Tbow shock,  an  expansion  fan  forms  at  point  g  and  the  freestream  flow  properties  are  used  as  input 
to  the  expansion  fan  flow  equations.  The  wedge  angle  in  this  case  is  only  a  function  of  angle  of  attack. 
If  a  =  0,  then  the  flow  properties  in  region  3  are  freestream,  that  is,  V3  «  Voc^ps  =  Poo,  and  a 3  =  a0 Q. 
If  a  >  0,  a  shock  forms  at  point  g  with  the  shock  angle  computed  using  a  wedge  angle  of  a.  If  a  <  0, 
an  expansion  fan  forms  on  the  underside  of  the  engine  and  flow  properties  in  region  3  are  calculated 
using  a  wedge  angle  of  —  a. 

For  the  control  surface,  the  flow  behind  the  leading  edge  of  the  elevator  is  determined  by  the  elevator 
deflection  angle  and  the  angle  of  attack.  More  specifically,  if  Se  =  —a  then  both  the  top  and  bottom  of 
the  elevator  experience  the  freestream.  Therefore,  V4  —  V5  «  Vqo,  Pa  —  p$  =  Poo,  and  —  <25  =  a^.  If 
Se  >  —a,  then  an  expansion  fan  forms  on  the  top  of  the  elevator,  while  the  bottom  of  the  elevator  experiences 
compression  and  a  shock  forms.  In  either  case,  the  wedge  angle  is  a  +  5e.  If  Se  <  —  a,  a  shock  forms  on  top 
of  elevator  and  an  expansion  fan  is  on  the  bottom  of  the  elevator.  In  this  case,  the  wedge  angle  for  the  shock 
and  expansion  fan  is  —  a  —  Se. 


VIII.  Total  Forces  and  Moments 


Having  determined  the  forces  on  each  of  the  surfaces,  the  moments  about  the  c.g.  that  each  force  produces 
must  be  determined.  To  do  this,  the  location  of  each  force  on  the  vehicle  must  be  computed.  Figure  3  shows 
the  forces  acting  on  the  vehicle.  Consider  first  the  upper  surface  with  the  forebody  and  aftbody  forces  given 


C 


f 


g 


h 


F, 


Igh 


mg 


Figure  3.  Forces  Acting  on  the  Hypersonic  Vehicle. 


in  Eq.  39.  The  forebody  force  acts  at  the  point 


(71) 
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while  the  upper  surface  afterbody  force  acts  at 

Fcfa  '■  (- y  >  -  (xf  +  y)  ta  nn,[/)  (72) 

For  the  lower  surfaces,  the  forces  act  at 

Fed  ■  (xf-  - y ,  y  tann,ij  (73) 

Fs/>  :  (xf  ~  Lf  ~  ^*i/tanriiL  +  hi)  (74) 

The  force  due  to  the  external  nozzle  acts  at  a  point  given  by 

F ef  ■  {xf  -Lf-Ln-  xFef  cos  (rltu  +  r2) ,  Lf  tan  njL  -  xFe/  sin  (ThU  +  r2))  (75) 

where  ^Fe/  is  the  x  point  of  the  center  of  mass  of  the  pressure  distribution  on  the  rear  ramp  assuming  an 
axis  system  centered  at  point  e  with  x  axis  pointing  along  line  ef  and  the  z  axis  pointing  up.  In  other  words, 
the  center  of  mass  of  the  pressure  distribution  was  computed  in  a  local  coordinate  frame  (local  to  surface 
ef).  Then,  this  distance  was  referenced  to  the  c.g.  of  the  vehicle.  In  Eq.  75,  x-pef  is  given  by 


.  _  L« 

•^Fe  f  ,  x  ,  p 

cos  (tijC7  +  t2)  In 

The  elevator  force  acts  at 

Fe  :  (  xcs ,  zcs )  (77) 

Equations  39,  40,  44,  63,  and  64  give  the  axial  and  normal  force  components  of  the  forces  acting  on  the 
vehicle.  Now,  the  moments  due  to  these  forces  can  be  calculated.  Positive  moment  is  defined  as  clockwise  or 
the  direction  that  tends  to  increase  angle  of  attack,  while  negative  moment  is  defined  as  counter-clockwise. 


The  moment  arms  are  given  by  Eqs.  71  -  75,  77. 

MCf,  =  (!  -  tan2  n,u)  (78) 

-  -2 

Mcfa  =  -Pixatan2nif/  (xf  +  y)  +Fiy  (79) 

Mcd  =  -P2Lf  tanri,L  (L/ ta^lri^  )  +  p2Lf  (y  _  ^)  (80) 

Mgh  =  P3Ln  (xf  -  Lf  -  ^  (81) 

Mef  —  ±Xef  [Lf  tanrijL  -  xFe/  sin  (t^u  +  r2)]  -  Zef  [xf  -  Lf  -  Ln  -  xFe/  cos  (t^u  +  r2)]  (82) 


where  xFe/  is  given  in  Eq.  76  and  Xef,Zef  are  expressed  in  Eq.  44.  The  sign  on  the  first  component  of  Mef 
will  depend  on  the  vehicle’s  geometry.  The  rule  used  to  determine  the  sign  is  as  follows: 


+  if  Lf  tanri^  —  xFe/  sin  (ritu  +  r2)  >  0 

-  if  Lf  tanri.L  -  xFef  sin  (ntu  +  t2)  <  0 

The  moments  produced  by  the  elevator  are 

Meu  =  -P$Le  sin  Sezcs  +  P5Le  cos  5excs  (84) 

MeL  =  P4Le  sin  Sezcs  -  P4Le  cos  5excs  (85) 

The  total  aerodynamic  forces  and  moments  on  the  vehicle  are 

Xtotai  =  Xcff  +  Xcfa  +  Xcd  +  Xgh  +  Xef  +  XeL  +  Xeu  (86) 

Z total  —  zcfs  +  Zcfa  +  Zcd  +  Zgh  +  Zef  +  ZeL  +  Zeu  (87) 

M total  =  Mcff  +  Mcfa  +  Mcd  +  Mgh  +  Mef  +  MGl  +  Meu  (88) 
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IX.  Stability  Derivatives 


The  stability  derivatives  of  interest  for  this  vehicle  configuration  are  Cza,  Cxot,  Cm  a,  Czq>  and  Cmq- 
These  computations  are  complex  and,  in  an  attempt  to  simplify  them,  the  following  notation  will  be  used 
when  necessary: 


L 

L 

L 

L 


_  (-Zvq  3?q) 

e/2  —  2 

_  (-^q  X  a ) 

e/s  —  3 

(Lg-Xg) 
e/4  4 

_  Xa ) 

e/5  —  5 


2 


3 


4 


5 


(89) 


.•I,  =  peae  +  (pe  Poo)(ae  A  OO  )(^a-^a)2 

^  a 

{pe  -  poo)  +  fy  (ae 

A3 


-t/q  ^°°)  ^a)  £a  (^e  ^00)  (-^a  ^a) 

^00)  ]^2  (yPe  Poo)  {p>e  ^00)  (I'a  *^a) 

(/>e  ^oo  )  (®e  ®oo) 

LI 


(90) 


fi  =  tan(ri>c/  +  r2) 
f‘2  =  itailTij; 

/3  =  sin  (ri,q  +  r2) 

U  =  cos  (ri,q  +  r2) 

Additionally,  the  pressure  on  the  rear  ramp  of  the  vehicle  (surface  ef)  is  given  by  Eq.  41,  which  is  a  function  of 
the  distance  moved  along  the  rear  ramp.  Thus,  this  pressure  distribution  is  not  constant  along  the  surface. 
At  point  e,  the  pressure  is  given  by  the  engine  exit  properties  (Pe),  while  at  point  f  the  pressure  is  the 
freestream  pressure  (Poo)-  An  obvious  choice  for  the  temperature  distribution  on  the  rear  ramp,  Tey,  is  a 
form  similar  to  the  pressure,  so  that 


Tef(s) 


1  + 


S(T1,[/  +  T2) 


(92) 


Then,  the  speed  of  sound  and  density  can  be  calculated  from  the  definition  of  speed  of  sound  and  the  perfect 
gas  law: 


aef  —  \j  T-PTe/ 


Pef  — 


RTef 


(93) 


Both  aef  and  pef  factor  into  the  stability  derivative  calculations.  Unfortunately,  the  expressions  in  Eq.  93, 
when  used  in  the  stability  derivative  calculations  do  not  allow  determination  of  a  closed-form  solution.  In 
order  to  facilitate  a  closed- form  solution,  the  following  approximations  for  the  speed  of  sound  and  density 
on  the  rear  ramp  will  be  used  in  the  stability  derivative  calculations: 


Pef  =  (P”LaPe)  (X  -  {La  -  Xa})  +  pe 
aef  =  (a!!£aae)  ( X  -  {La  -  5a})  +  ae 


(94) 


where 

%a  —  %  —  L/a  Xa  (95) 

These  are  first-order  approximations,  which  capture  the  boundary  conditions. 


A.  a  Derivative  of  Z-Force  Coefficient 

To  compute  the  change  in  Z-force  component  due  to  a  change  in  angle  of  attack,  the  infinitesimal  force 
expression  that  contains  the  vertical  velocity  perturbation,  w,  must  be  integrated.  Thus, 


(cz), 


1 

Qoo$ 


rx  f  L  f  pLa—xa 

/  (dFgh)z_w  +  /  (dFef)z 

Jxf—Lf—Ln  J  —Xn. 


(96) 
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where  f  (dFcf)  _  is  the  differential  force,  on  the  upper  surface,  in  the  z  direction,  due  to  w  motion  and 
likewise  for  the  other  surfaces.  Using  the  appropriate  differential  force  elements  from  Eq.  25  in  Eq.  96 
produces 


+  - 


(Cz), 

1 

QooS  _ 

1 

PXf~Lf 

^00-5 

Jxf-Lf-L 7 

nxf  nxf 

/  —  piCliW  COS  Ti^udx  —  /  p2d2W  COSTi^dx 

j  —xa  Jxf—Lf 


-p3d3wdx  + 


L 


—  pefdefW  COS  JJ  +  T2)  dx 


Performing  the  integrations  yields 


1 


(Cz)w  =  [~Plal  COSTi tU  (xf  +  Xa)  ~  p2d2  COS  TliLLf  ~  p3d3Ln] 

-weos(ruu  +  r,)  +  +  AM 


Qoo  S 


If  <  1  then  y~  «  a.  Thus,  w  =  V^a.  Therefore, 


dCz  Ux 


da  QooS 


[~Pldi  COS  Ti^zjL  -  p2d2  COST\iLLf  -  p3d3Ln] 


Vqo  COS  (jxu  +t2) 
QooS 


[A\La  +  A2Lef2  +  A3L 


ef3\ 


B.  a  Derivative  of  X- Force  Coefficient 

The  change  in  X-force  due  to  a  change  in  angle  of  attack  can  be  calculated  using 

(cx)w = -E  [  0 dF)x_w = -E  n  (dFcf)x_w  +  r  (, dFcd)a 

Qoo>3  J  (loo^  J—  xa  Jxf—Lf 

pXf—Lf  pLa—xa 

/  (dFgh)x_w  +  /  (dFe/)a 

Jxf—Lf—Ln  J —Xn 


1 

QoqS 


Substituting  in  the  differential  forces  yields 

1 


+ 


(Cx)w  = 

1 


QooS 


Qoo  S 

xf~Lf 


rxf  rxf 

/  pidiw  sinri^jjdx  —  /  p2a2re  sinri^dx 

J  3?q,  J  Xf  Lf 


pX  j  Li  j  pLa—Xa 

/  —  p3d3wOdx  +  /  PefdefW  sin  (ti?£/  +  T2)  dx 

JXf—Lf—Ln  j —Xn 


,Xf~Lf~L 


(97) 


(98) 


(99) 


(100) 


(101) 


Performing  the  integrations,  simplifying,  and  assuming  that  y-  C  1  so  that  y-  «  a  gives 

1  W  r  •  r  •  r  1  ,  Voa  sin  (t^u  +t2)  (  . 

,Voo  [pidi  smr^f/L  -  p2a2  sinn^L/J  + - ^ - [^iLa  +  A2Le/2  +  A3Lef3 \  (102) 


da  q^S 


QooS 


C.  a  Derivative  of  Pitching  Moment  Coefficient 

For  Cm  a:  first  find  the  contribution  to  the  pitching  moment  due  to  a  velocity  w: 
1 


(C.\/),„  = 


qooSc 

1 

q^oSc' 

1 

qooSc2 


J  z(dF)x_w-J  x{dF)z_ 

/Xf  pXf  pXf—Lf  pjja—xa 

{dFcf)x_w  +  /  (dFcd)x_w  +  /  (dFgh)x_w  +  /  (dFe/) 

-CCa  Jxf—Lf  Jxf—Lf  —Ln  J—  xa 


(*L  a  Xa 


/xf  rxf  rxf~Lf  r 

(' d^cf)z_w  +  /  (' dFcd)z_w  +  /  (^Fsft)z_UJ  +  / 

-*a  Jxf—Lf  Jxf—Lf—Ln  J  — 


nLa—xa 


(dFef)z 


(103) 
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Substituting  the  appropriate  expressions,  performing  the  integrations,  simplifying,  and  assuming  that  y-  <C 
1  so  that  y-  «  a  yields 


1  1 


-hoc  [-pia  1  tanri,!/  sin  ti^L2  +  p2a2  tanri?L  sin 


dCM 

da  q00Sc2 

+  °°  — ^  A- - ~  [/i  (AiLef2  +  A2Lef3  +  AzLef4)  +  (fixa  —  /2)  (-4iLa  +  A2Lef2  +  ^h^e/s)] 

qooDC 

H - [piaicosri,!/  (z/  -  £„)  +  p2£i2  cos  rltLLf  (2xf  -  Lf)  -  p3a3Ln  (Ln  -2xf  +  2Lf)\ 

q^oc  z 

,  KoCOS  (ThU+T2)  lA  T  ,  4  r  ,  4  r  ! 

■+  gg  [^l-t'e/2  +  MLef3  +  ^3-^e/4J 


(104) 


D.  q  Derivative  of  Z-Force  Coefficient 

The  Z-force  coefficient  due  to  a  pitching  motion  is 


<^ = ^  /  i  [/i  + c, 

pXf-Lf  rLa- 

/  0 dFgh)z_q  + 

J  Xf  —Lf  —  Ln  J  —  xa 


(dFcd)z_ 


z-q 


(dFef)z_q 


Performing  the  integrations  and  simplifying  yields 


+ 


Qoo  s 


dCz  1  L  4  . 

[cos  t\  jj  -  tann^sinn^L 

dq  q^oS  2 

+  ~~gP2a2~2  [C0STl>£  (2S/  ~  Lf )  —  tan ri,L  sin 'Ti.i-fi/] 

+~^~sf}3a3~^L  2S/  _  2-^/] 

[(/1/3  —  fi)  (AiLef2  +  A-2Lef,,  +  AzLef4 )  +  ( fif3Xa  —  f2fz)  {A\La  +  A2Lef2  +  A2,Lef3)] 


(105) 


(106) 


E.  q  Derivative  of  Pitching  Moment  Coefficient 

The  pitching  moment  due  to  a  pitch  rate  can  be  calculated  using 


z-q 


1 


qooSc 

1 

q^Sc 


(CM)q  =  J  Z  (dF)x_q  -  J  X  (dF) 

/x  f  pXf  pXf  L  f  pLa  xa 

(dFcf)x_q+  _  (dFcd)x_q+  (dFgh)x_q  +  (dFef)x_q 

-xa  J  x  f  Lf  J  x  f  Lf  Ln  J  xa 

/Xf  pXf  pXf—Lf  pLa-xa 

(dFcf)z_q  +  /  ( dFcd)z_q  +  /  ( dFgh)z_q  +  /  ( dFef)z_q 

-Xa  JXf—Lf  Jxf— Lf—Ln  J  — 


'Xf-Lf 


'  -Xa 


(107) 
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Substituting  the  appropriate  expressions,  performing  the  integrations,  and  simplifying  produces 


d  a 


M 


dq 


qooSc 


pidi  tarn  T\  u 


L 3  cosTir/  /  _o  ^_o 

tanri  (/  sinri  [/  —  -\ - ‘  ~  1 

3  6 


(~x3f  +  2xl  +  3  x2aXf) 


- 7^/9202  tan  Ti  l 

qooSc 

fi 


L3f  Lf 

tanri)L  sinri)L-^-  -  cost^t,-^-  (2  Xf  -  Lf) 


-  [( A\Lef3  +  A2Lef4  +  A3Lef5)  (fif3  +  /1/4)] 


fi 


QooSc 


QooSc 

\(AiLef2  +  A2Lef3  +  A^Lgf^)  (2/j  fzxa  —  2/1/2/3  —  f2fA  +  fif4%a)] 


fi 


QooSc 


-  \{AiLa  +  A2Lef2  +  AsLef3 )  (fif3Xa  —  2flf2f3xa  +  /I/3)]  (108) 


qooSc 


~P  i«i 


tan  T\  u  sin  T\  ij 


(-Xf  +  3 x2aXf  +  2a;„)  H - (x3f  +  x®) 


qooSc 


~P  2«2 


tanri^  sinri^ 

6 


L2f  (2 Lf  -  3xf)  +  cos 7~i5/, ~ —  (3 x2f  -  3 XfLf  +  L2) 

O 


q^Sc 


( xf  Lf)  ( xf  Lf  Ln) 


qooSc 


-  [( A\Lefz  +  A2LefA  +  A3Lef5)  (/1/3  +  /4J  +  ( A\Lef2  +  A2Lef3  +  A3Lef4 )  ( fifsxa  —  f2fs)} 


X.  Engine 


For  hypersonic  flight,  the  propulsion  system  necessary  to  produce  the  required  thrust  is  either  a  rocket  or 
a  supersonic  combustion  ramjet  (scramjet).  The  advantage  of  a  scramjet  over  a  rocket  is  that  the  scramjet 
is  an  airbreathing  propulsion  system,  thus  eliminating  the  need  to  carry  the  oxidizer  onboard  the  vehicle. 
This,  in  turn,  allows  for  increased  payload.  The  scramjet  model  used  in  this  work  is  identical  to  that  used  by 
Chavez  and  Schmidt.5  Engine  inlet  conditions  are  primarily  determined  by  the  flow  behind  the  oblique  shock. 
The  scramjet  consists  of  3  sections;  a  diffuser,  a  combustor,  and  an  internal  nozzle.  The  flow  through  the 
diffuser  and  nozzle  is  assumed  to  be  isentropic,  quasi-one-dimensional,  while  the  flow  through  the  combustor 
is  assumed  to  be  quasi-one-dimensional  in  a  constant  area  duct  with  heat  addition.  The  working  fluid  in  the 
engine  is  assumed  to  be  a  perfect  gas  with  constant  specific  heats.  There  are  two  control  variables  which 
affect  the  engine:  diffuser  area  ratio,  Ad,  and  temperature  addition  in  the  combustor,  ATo.  Figure  4  shows 
the  engine  model. 

Inlet  conditions  to  the  diffuser,  which  are  inlet  conditions  to  the  engine  module,  are  determined  from 
the  flow  analysis  of  the  vehicle’s  lower  forebody  and  vehicle  geometry.  The  lower  forebody  flow  is  turned 
parallel  to  surface  cd.  To  determine  the  engine  inlet  conditions,  the  flow  must  be  turned  parallel  to  the 
engine,  with  the  turning  angle  given  by  Oblique  shock  relations  are  used  to  determine  the  engine  inlet 
flow  properties,  with  the  flow  properties  in  region  2  being  the  input,  and  Me.n,Pe.n,  and  Te.n  (the  engine 
inlet  flow  properties)  being  the  outputs  of  the  flow  calculations.  By  turning  the  flow  to  be  parallel  to  the 
engine,  a  force  and  moment  is  imparted  on  the  vehicle.  The  force  is  given  by 


Fxinlet  =  lM$P2  (1  -  cos  (ri)L  +  a)) 
Fernet  =  lMiP2  sin  (n,L  +  a) 


This  force  acts  at  point  d  in  Figure  1.  Hence,  the  force  acts  at 


^inlet 


Xf  —  Lf  L/tanri?£ 


Thus,  the  moment  produced  is  given  by 


Miniet  —  Lf  ta,nTi^LFXirilet  (xf  Lf)FZinlet 


(109) 


(110) 


(111) 
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Diffuser 

Combustor 

^  fcJ 

Internal  Nozzle 

^  w 

^  w 

^  w 

Figure  4.  Scramjet  Engine. 


Now,  the  flow  must  be  propagated  through  each  section  of  the  scramjet  to  determine  the  engine  out¬ 
put  properties.  For  the  diffuser,  continuity  is  applied  to  determine  the  Mach  number  at  the  diffuser 
exit /combustor  inlet: 


[1  +  1(7— l)M2](^)  [1  +  1(7  — i)M2j(^) 

Ml  D  Ml. 


(112) 


where  Mc  is  the  Mach  number  at  the  combustor  inlet  and  Mein  is  the  Mach  number  at  the  engine  inlet.  The 
pressures  and  temperatures  at  the  combustor  inlet  are  given  by 


Pr  =  P, 


Tr=Te 


l  +  §(7- l)Me2in 
1  +  A(7-1)Mc2 

l  +  A(7-l)M2n 
1  + |(7-1  )M2 


(113) 


(114) 


where  Pc  and  Tc  are  the  pressure  and  temperature  at  the  combustor  inlet,  and  Pein  and  Tein  are  the  pressure 
and  temperature  at  the  engine  inlet.  For  the  combustor,  the  exit  Mach  number,  temperature,  and  pressure 
are  calculated  using 


Ml[  1  +  |  (7  -  1  )M2]  _  Ml  [1  +  A  (7  -  1)MC2] 


(7  Ml  +  1)J 


P  —  P 

1  n  —  J-  c 


(7  Ml  +  l)2 

7  Ml  +  1 

7M2  +  1 


+ 


Ml 


ATn 


(7Mc2  +  1)' 


=  TflK± 


1  M„ 


V7  Ml  +  1  Mc 


(115) 

(116) 
(117) 
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where  Mn,  Pn,  and  Tn  are  the  Mach  number,  pressure,  and  temperature  at  the  inlet  to  the  nozzle,  respectively 
and  ATo  is  the  increase  in  total  temperature  across  the  combustor  due  to  the  combustion  of  fuel.  For  the 
nozzle,  the  exit  properties  are 


[1  +  |(7  — 1  )Me2](?-)  _  [1  +1(7-1)M„2]^) 

Ml  N  Ml 


(118) 


where  Me  is  the  Mach  number  at  the  engine  exit.  The  pressures  and  temperatures  at  the  engine  exit  are 
given  by 


P  —  P 

±  e.  -1  n 


T,  =  T, 


1+  i(7-l)M2 

1.1  +  5  (7  —  1)  Ml 

1+ 1(7-1  )Ml 


Ll  +  1(7-1)M|J 


(119) 

(120) 


where  Me,  Pe,  and  Te  are  the  Mach  number,  pressure,  and  temperature  at  the  engine  exit,  respectively,  and 
An  is  the  internal  nozzle  ratio  defined  as  the  ratio  of  the  nozzle  exit  area  to  the  nozzle  inlet  area. 

The  thrust  per  unit  width  of  the  engine  module  is  given  by4 


T  —  fn0  {Ve  Foo)  +  (pe  Poo) 


Ae 

b 


( P2  “Poo) 


Ad  An 


(121) 


where  T  is  the  engine  thrust,  m0  is  the  mass  flow  through  the  engine,  Ve,  pe  are  the  flow  velocity  and  pressure 
at  the  engine  exit,  V25P2  are  the  flow  velocity  and  pressure  at  the  engine  inlet,  is  the  freestream  pressure, 
Ae/b  is  the  exit  area  per  unit  width,  Ad  is  the  diffuser  area  ratio,  and  An  is  the  nozzle  area  ratio.  The 
moment  produced  by  the  thrust  force  is 


M-engine  —  ^PjtanTi^  ^  )  T 


(122) 


The  mass  flow  through  the  engine  is  a  function  of  the  shock  angle.  Essentially,  in  order  to  calculate  ra0, 
it  is  necessary  to  calculate  how  much  mass  flow  the  engine  captures.  Figure  5  sets  up  the  geometry  for  this 
calculation.  It  is  assumed  that  the  vehicle  is  operating  such  that  a  shock  forms  off  the  forebody  surface.  Let 
the  engine  inlet  capture  area  be  denoted  by  Aq  and  the  spill  area  be  denoted  by  As.  When  a  shock  forms, 
use  Figure  5  to  compute  the  capture  area,  which  becomes 


A  _  (L.f-HLf-(L.f  ^Ti^L+hi)  cot ((3icd- a)})  p 
^0  ~  cos (Plcd-a)  Sm  Pled 

where  /3icd  is  the  shock  angle  for  lower  surface  cd  and  9\  is  defined  as 

_i  f  Lf  tan ti,l  +  hi 


0 1  —  tan 

The  mass  flow  through  the  engine  becomes 


7Ti0  —  Pqo  TTq, 


RTrs, 


-An 


(123) 


(124) 


(125) 


With  the  inclusion  of  the  stability  derivatives,  the  engine  inlet  turning  force  and  moment,  and  the  thrust 
and  resulting  moment,  the  total  aerodynamic  forces  and  moments  on  the  vehicle  are 


Xtotai  —  Xcff  +  Xcfa  +  Xcd  +  Xgh  +  Xef  +  XeL  +  X&u  +  q^S  ^  a.  +  FXinlet 


(126) 


Ztotal  —  Zcff  +  Zcfa  +  Zed  +  Zgh  +  Zef  +  ZeL  +  Zejj  +  QooS  ^  Q  +  QooS  - f  FZinlrf  (127) 


dq  2W 


Mtotai  —  Mcff  T  Mcfa  Med  T  Mgh  -\-  Mef  -\-  MeL  -\-  Meu  -\-  qooSc  a  T  qooSc  T  Miniet  -b  Mengine  (128) 

J  da  dq  2  Voo 
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Figure  5.  Geometry  for  Mass  Flow  Capture  Area. 


XI.  Aeroelastic  Effects 

Thus  far,  only  a  rigid  body  has  been  considered.  In  fact,  the  vehicle  is  long  and  slender  and  thus  highly 
flexible.  This  structural  bending  affects  downstream  flow  and  should  be  incorporated  into  the  model.  In 
this  section,  the  aeroelastic  effects  are  included  in  the  analysis.  Piston  theory  is  still  used  to  determine  the 
pressure  distribution  on  the  surfaces  of  the  vehicle  and  much  of  the  analysis  already  presented  can  be  easily 
adapted  to  included  these  additional  effects. 


A.  Aeroelastic  Model 

In  order  to  develop  the  aeroelastic  model,  a  few  assumptions  are  made.  First,  the  flexible  vehicle  is  modelled 
as  a  cantilever  beam  fixed  at  the  c.g.  (only  the  forward  section  of  the  vehicle  is  currently  considered  flexible). 
Second,  the  beams  are  assumed  to  have  constant  mass  density,  area,  and  flexural  rigidity  (El),  where  El  is 
chosen  to  give  the  desired  natural  frequency  of  vibration.  Also,  it  is  assumed  that  the  flexible  effects  only 
perturb  the  surface  velocities  in  the  z  (normal)  direction.  This  assumption  is  justified  using  the  small  angle 
approximation,  i.e.,  the  deflection  of  the  tip  of  the  beam  is  small  compared  to  the  length  of  the  beam.  Lastly, 
it  is  assumed  that  the  pressure  distribution  over  the  forebody  of  the  vehicle  is  constant  even  when  the  beam 
is  deflected. 

The  transverse  vibrations  in  the  beam  satisfy  the  following  partial  differential  equation:6 

d4w(x,t)  ^d2w(x,t) 

where  w(x,t)  describes  the  position  of  the  beam,  E  is  Young’s  Modulus,  and  I  is  the  moment  of  inertia  of 
the  beam  cross-section  about  the  y-axis.  To  solve  this,  use  the  method  of  separation  of  variables.  Assume 


(129) 


w(x,  t )  =  $(x)T(t) 

Substituting  the  expression  for  w(x,t)  in  Eq.  130  into  Eq.  129  and  simplifying  yields 

El  <94$(x)  _  1  d2T(t) 

m$(x)  dx 4  T(t)  dt 2 


(130) 


(131) 
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Since  the  left  side  of  Eq.  131  does  not  change  as  time  varies,  the  right  side  of  Eq.  131  must  be  a  constant. 
Similarly,  since  the  right  side  of  Eq.  131  does  not  change  as  x  varies,  the  left  side  of  Eq.  131  must  be  a 
constant.  Let  this  constant  be  cj2,  such  that 


El  dA$(x)  _ _ 1  d 2T{t)  _  2 

ra<L(x)  dx 4  T(t)  dt 2 


(132) 


Utilizing  the  method  of  separation  of  variables,  Eq.  132  can  be  written  as  two  differential  equations,  one 
with  respect  to  position  and  one  with  respect  to  time: 


<94<f>(x) 

dxA 


/34$(£)  =  0 


where  /34 


d2T{t) 
dt 2 


+  uj2T(t)  =  0 


The  general  solution  to  Eq.  133  is6 

$(#)  =  Ci  sin  (3x  +  C2  cos  kx  +  C3  sinh  kx  +  C4  cosh  kx 


(133) 

(134) 


(135) 


Using  the  following  boundary  conditions  for  the  forward  beam 


d2${x)  I 

dx  2  \x  =  X 


$(®)U=0 


d3$(x) 

dx3 

d$>(x) 

dx 


f  —  0 
U=o  =  0 


(136) 


which  state  that  the  bending  moment  and  shear  force  are  zero  at  the  free  location  (x  —  Xf)  and  the 
displacement  and  slope  are  zero  at  the  fixed  location  (x  —  0),  along  with  the  modal  shape  expression, 
(Eq.  135),  and  simplifying  results  in  the  frequency  equation 


cos  /3L  cosh  /3L  =  —  1 


(137) 


Eq.  137  has  an  infinite  number  of  solutions  with  the  first  few  given  by 

pL  =  1.8751, 4.6941, 7.8548, 10.9955, 14.1372, . . . 


(138) 


The  values  of  /3r,  r  =  1,2,3, .. .  in  Eq.  138  are  called  the  eigenvalues.  Corresponding  to  these  eigenvalues, 
the  natural  modes  of  the  forward  beam  are6 


$f^r(x)  =  Af)V  [(sin /3f,rXf  —  sinh /3f,rXf)  (sin (3f,rx  —  sinh /3f,rx)] 
-\-Af^r  [(cos  (3f,r%f  +  cosh  /3f,rXf)  (cos  /3fjVx  —  cosh  f3fjrx)] 


where  Af^r  is  a  normalizing  factor,  selected  such  that 

rxf 


rxf 

J  m^2r{x)dx  M  1 


Thus,  Af,r  becomes 


-T/,r  — 


m  A 


fxp,  +  A/>P2  +  ^4/,r  p3 


where 


(sin  f3fcxf  —  sinh  (3  fc  x 


4  Pk 


*f,rp2 


[-2  cos  pkxk  sin  f3kxk  +  sinh2/3fca:fe  -  4sin/3kxk  cosh  f3kxk  +  4  cos  f3kxk  sinh/3fea;fc] 


(sin  (3k  x  j  —  sinh  /3kx  [cos  /3k  x f  +  cosh  )3kx ^ 


Pk 


(sin  (3kxf  -  sinh/3kxfy 


(cos  PkXf  +  cosh  Pkx  f')2 


rP3 


4  Pk 


[2  cos  (3kxk  sin  /3kxk  +  sinh  2/3kxk  —  4  cos  /3kxk  sinh  Pkxk  —  4  sin  (3kxk  cosh  /3kxk  +  4(3kxk] 


(139) 

(140) 

(141) 

(142) 

(143) 

(144) 
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B.  Forced  Response 


Let  the  forcing  function  in  Eq.  129  consist  of  distributed  and  concentrated  loads  so  that  Eq.  129  can  be 
written  as 


El 


d4w(x,  t )  ^  d2w(x ,  t) 


+  m- 


dx 4  dt 2 

where  5(#)  is  the  dirac  delta  function  defined  as 


=  f(x,t)  +  Fj(t)S(x  -  Xj) 


(145) 


V!<x7‘ 

{  0  if  x  ±  xj 

From  the  expansion  theorem,  the  solution  to  Eq.  145  is 

oo 

w{x,t)  = 

r= 1 


(146) 


(147) 


where  rjf,r(t)  is  the  generalized  modal  coordinate,  for  the  forebody  beam,  that  satisfies 

Vf,r  (t)  +  2C  f,rUf,rVf,r  (*)  +  ^2f,rVf,r  (t)  =  Nf,r  ( t )  (148) 

Here,  is  a  generalized  force  for  the  kth  mode  shape  of  the  forebody  beam,  defined  by6 


rx  f  n 

Nf,r(t)  =  /  ®f,r(x)f(x,t)dx  +  y^^f,r{Xj)Fj{t) 

J0  ,=1 


(149) 


where  n  is  the  number  of  concentrated  loads  on  the  beam.  Given  the  loading  on  the  forebody  beam,  the 
generalized  force  for  the  first  mode  becomes 


NfA(t)=  f  $ftl(x)Pcfdx-  f  $fti(x)Pcddx 
JO  J  Xf—Lf 


(150) 


Evaluating  this  expression  produces 


Nf,i(t)  -  LiA-Sl  [(—  sin  {0iXf}  +  sinh {/?!*/})  (cos{/?ia;/}  +  cosh{/3i®/}  -  2)] 


Pi 


Af,ipcf 

Pi 


[(cos{/?ia;/}  +  cosh  {PiXf})  (sin  {PiXf}  -  sinh  {PiXf})] 


AjjP 


cd 


Pi 

Af^\Pcd 


[(-sin{/3i  Xf}  +  sinh{/3ix/})  (cos  {0iXf}  -  cos  {/3i  (xf  -  Lf )}  +  cosh  {/3i  Xf}  -  cosh{/3i  (xf  -  Lf)})] 

[(cos  {/?i  Xf}  +  cosh  {/?i  Xf})  (sin  {PiXf}  -  sin  {f3i  (xf  -  Lf)}  -  sinh{/3i  Xf}  +  sinh{/3i  (xf  -  Lf)})] 

(151) 


The  solution  to  Eq.  148  for  the  rth  mode  is6 

1 

Vf, r(t)  — -  /  Nftr(r)  exp_c-f'r“-f'r(t_T)  sin  (wftrd(t  -  r))  dr+ 

^f,rd  J  0 


exp 


-C  f,rUf,rt 


COS  UJf^rdt  + 


C/,i 


—  sin  ujf^rdt 


VfA  0)  + 


^ f,rd 


exp 


-C  f,rVf,rt 


(X  tf,r) 

where  Wf,rd  —  w f,r  (l  —  C/r)  ■  The  initial  conditions  are  given  by6 

77/>r.(0)  =  Jq*  m&f:r(x)w(x,  0)dx 
VfA0)  =  f0Xf  r(x)w(x,  0 )dx 


sin  u)firdt 


(152) 


vfA  o) 


(153) 
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Evaluating  the  expressions  in  Eq.  153,  for  the  first  forebody  mode,  yields 


(154) 

In  order  to  incorporate  aeroelastic  effects  into  the  model,  a  few  simplifying  assumptions  are  made.  First, 
the  vehicle  does  not  stretch  or  compress  along  the  x-axis.  Second,  for  small  displacements,  when  the  beams 
(vehicle)  flexes,  there  is  no  change  in  the  x  direction  displacement.  With  these  assumptions,  aeroelastic 
effects  only  occur  in  the  z-direction. 

The  aeroelastic  effects  can  be  accounted  for  by  taking  the  time  derivative  of  Eq.  147, 

oo 

w(x,  t)  =  (155) 

r= 1 

and  including  this  effect  in  the  expressions  for  the  velocities  on  the  upper  and  lower  surfaces,  namely  Eqs.  11, 
14,  and  15.  The  differential  forces  on  the  upper  and  lower  forebody  surfaces,  taking  into  account  the  forward 
section,  becomes 


dFcf  =  (—Pi  —  Picli  {[u  +  q  tanri?£/  ( x  —  Xf)\  sinri?t/  —  [w  —  qx  +  w(x,  t )]  cosri^})  dAcfncf  (156) 

dFcd  —  (-P2  -  p2a2  {[ u  -  gtanri?L  (x  -  xf)]  sinr1?L  +  [w  -  qx  +  w(x,  t)]  costal})  dAcd ncd  (157) 

With  these  differential  forces,  stability  derivatives  due  to  bending  of  the  vehicle  can  be  determined.  For  the 
normal  force, 


(Cz)^  = 


Qoo^ 


rxf  rxf  1  rxf  rxf 

/  (dFcf)z_ih+  {dFcd)z_lj]  = - -  -p1a1wcosT1,udx+  -p2a2w  cos  r^Ldx 

JO  J Xf—Lf  Qoo ^  Jo  Jxf—Lf 

(158) 


Substituting  Eq.  155  into  Eq.  158  produces 


(CZ)yj  ~ 


Qoo^ 


rxf  rxf 

/  (dFcf)z_th+  _  (, dFcd)z_v 

JO  J  Xf  —Lf 


I  rx  f  00  cx  f  00 

= - n\  -Piaiy2$f,r{x)r)f,r{t)cosTitUdx  +  ~P2d2  T 

do  r=l  Jxf-Lf  r=1 


(159) 


(x)riftr(t)  cosTi^dx 


which,  for  one  mode  becomes 
dCz 


z  1  rxf  rxf 

0.  —  =  - 77  /  -piai4>/,i(x)  cosri, udx  +  /  -p2a2®  f, i(x)  cos  n,Ldx 

or]f,i  qooO  Jo  Jxf—Lf 

(sm/3iXf  —  sinh/3iXf) 


Qoo  S 


-2piai  cosn,i/A/,i 


p2a2  cos  Ti^Af^i 
PiqooS 

p2a2  costi,l^4/,i 


(160) 


[(sm/3iXf  —  smh/3iXf)  (cos/3i  ( Xf  —  Lf)  +  cosh/^i  ( Xf  —  Lf))] 


Pi  qooS 


[(cos^i Xf  +  cosh^i Xf)  (—  sin^i  (xf  —  Lf)  +  sinh/3i  (xf  —  Lf))] 


where  Af^  is  the  normalizing  factor  associated  with  the  first  mode.  Using  the  assumption  that  the  flexible 
effects  only  perturb  the  surface  velocities  in  the  z-direction,  the  axial  force  stability  derivative  associated 
with  the  flexible  effect  becomes 

dCx  "  (161) 


<%,  1 


=  0 


For  the  pitching  moment, 


(CWU  =  ~ 


q^Sc 


rxf  rxf 

/  x  ( dFcf)z_dj  +  x  (dFcd)z_ti. 

JO  Jxf—Lf 


(162) 
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For  the  first  mode  only, 


dCm  _  ^Af  ipiai  cosTi’jj 


dlif,  i 
,1  p: 

P\ 

AfjiP2CL2  COSTi^l 


PlqooSc 


[l  —  cos  pi  Xf  +  cos  Pi  Xf  cosh  P\Xf  —  cosh  PiXf] 


H - f'lpft2  2  C-  1L  [2  -  cosPiLf  +  Pi  ( Xf  -  Lf)  sin  PiLf  +2  cosPiXf  coshpiXf] 

Pi  Qoo  L  c 


pfqooSc 


[pi  ( Xf  —  Lf )  ( sinpiXf  cosh^i  ( Xf  —  Lf)  +  cosPiXf  sinh/3i  ( Xf  —  Lf ))] 


H - f,1%  2  c-  1L  AnpiXf  sinh/3i  ( Xf  -  Lf)  -  cosPiXf  cosh^i  ( Xf  -  Lf)]  (163) 

P  i  q<x>Lc 

4 - ^  ^Pq2  2  g-  1L  [AnhpiXf  sin^i  ( Xf  —  Lf)  —  cosh/3i^/  cos/5i  ( Xf  —  Lf)] 

Pi  qoobc 

+  f'lPo22  C^Tl,L  [-pi  ( Xf  -  Lf)  sinh/3i^/  cos^i  ( xf  -  Lf)  -  Pi  (xf  -  Lf)  cosh/?ix/  sin^i  ( xf  -  Lf)] 

Pi  qoo&c 

+  yl/’ip2a2CosTi’L  ^  sinh/3lLf  _  cosh/Sii/] 

With  the  inclusion  of  the  flexible  stability  derivatives,  the  total  aerodynamic  forces  and  moments  on  the 
vehicle  are 

rfC1 

Xtotai  =  Xcff  +  Xcfa  +  Xcd  +  Xgh  +  Xef  +  XeL  +  Xeu  +  q^S  -^-cx  +  FXirilet  (164) 

Ztotal  =  Zcff  +  zcfa  +  Zed  +  Zgh  +  Zc/  +  ZeL  +  Z6[/  +  q^S-^-a  +  +  ^inZet  +  3ooS^T^?7/,l 


Mtotal  —  Mcf.  -\-Mcfa  +Mcd+Mgh+Mef+MeL  -\-Meu  -\-qooSc—~ - OL-\-qooSc—z - - VMinietp-Menginep-qooSc— - f]fp 

ocx  Oq  zVo o  Or]fp 

(166) 


XII.  Equations  of  Motion  for  a  Flexible  Vehicle 

The  equations  of  motion  for  the  flexible  vehicle  are4 

Vt  =  ^  (T  cos  a  —  D)  —  g  sin  (6  —  a) 

d  (“Tsin«  -  L)  +  Q  +  Vt  COS  ~  a) 

lyyQ  ~  ^ fVf  =  M 

kfVf  +  2  CfUffif  +  u2fVf 

h  =  Vt  sin  (6  —  a) 

where 

=  fof  x$fAx)dx 

kf  =  1  +  p- 

J  lyy 

Clearly,  the  flexible  effects  are  coupled  into  the  pitch  rate  equation.  In  addition  to  this,  the  bending  of  the 
structure  has  an  effect  on  the  angle  of  attack  of  the  vehicle.  Since  engine  performance  is  a  function  of  shock 
angle  and  shock  angle  is  a  function  of  angle  of  attack,  a  significant  change  to  the  vehicle’s  performance  can 
occur  due  to  structural  bending.  It  is  assumed  that  the  entire  forebody  observes  the  same  change  in  angle 
of  attack  as  seen  at  the  nose  of  the  vehicle.  In  other  words,  the  worst  case  change  in  angle  of  attack  is  used 
for  the  entire  forebody.  This  change  in  angle  of  attack  is  computed  as: 

A  a  =  arctan  [&  (xf)  rjf,i{t)]  (169) 

XIII.  Results 

At  this  point,  only  initial  simulation  results  are  available.  Currently,  the  model  has  been  simulated  open- 
loop,  with  a  fixed  control  surface  deflection,  to  ensure  that  the  model  is  operating  correctly.  One  point  of 
interest  that  can  be  obtained  from  this  simple  simulation  is  the  contribution  to  the  forces  and  moments  due 
to  the  inclusion  of  unsteady  effects.  Figure  6  shows  the  steady  and  unsteady  X  force,  Figure  7  shows  the 
steady  and  unsteady  Z  force,  and  Figure  8  shows  the  steady  and  unsteady  pitching  moment.  Obviously,  the 
unsteady  components  will  have  an  impact  on  the  total  forces  and  moments,  as  these  terms  are  not  negligible 


(167) 


(168) 
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compared  to  the  steady  terms.  In  fact,  the  unsteady  pitch  moment  has  the  largest  effect,  followed  by  the 
unsteady  Z  force.  The  unsteady  X  force  has  little  effect  on  the  total  X  force. 
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Figure  6.  Steady  and  Unsteady  X  Forces. 


XIV.  Conclusions 

In  this  work,  piston  theory  is  used  to  develop  a  model  for  the  longitudinal  dynamics  of  a  hypersonic 
vehicle.  In  particular,  velocities  of  flow  normal  to  the  surface  of  the  vehicle  are  used  in  a  first  order  piston 
theory  framework  to  determine  the  pressures  on  the  surfaces  of  the  vehicle.  The  pressures  are  then  integrated 
over  the  body  to  determine  the  forces  acting  on  the  vehicle.  Piston  theory  is  useful  here  because  it  allows 
the  inclusion  of  unsteady  aerodynamic  effects,  which  are  not  captured  using  other  techniques. 
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Figure  7.  Steady  and  Unsteady  Z  Forces. 

x  -fo4  Steady  and  Unsteady  Pitching  Moment  vs.  Time 


Figure  8.  Steady  and  Unsteady  Pitching  Moments. 
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